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Abstract 

We report the observation of the decay I?si(2536)"'" D~^'k~K~^ . We also measure the hehcity 
angle distributions in the decay -Dsi(2536) D*^Kg and thus constrain the contributions and the 
phase difference of D and S wave amplitudes in this decay. The results are based on a 281 fb~^ 
data sample collected with the Belle detector near the T(45) resonance, at the KEKB asymmetric 
energy e'^e" collider. 

PACS numbers: 
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INTRODUCTION 



Two states Dsj{2317)^ and Dsj{2A60)^ have been discovered recently both in continuum 
e'^e~ annihilation near -y/i = 10.6 GeV/c^ and in B meson decays 0- Their decay 
properties are consistent with the assumption that these are = 0^, 1^ states with j = 
L + Ss = 1/2. Here L = 1 is the orbital momentum, Sg is the spin of the light antiquark. 
However, their masses are unexpectedly low This has renewed interest in measurements 
of P-wave excited charm mesons. 

We report the first observation of the decay Dsi(2536)+ D^n^K^ (the inclusion of 
charge conjugate modes is implied throughout the paper). The D^n^ pair in the final state is 
the only Dir combination that cannot come from a D* resonance. Note that D*^ mesons can 
only be produced virtually here since M^j.o < M£)++M^-. The Dsi{253Q)~^ — >■ D^tt~K^ and 
^51(2536)+ Dfn^7r~ p modes are the only known three-body decays of the Dsi(2536)"'". 

In addition, we have performed an angular analysis of the Dsi{253Q)~^ D*^K^ mode. 
In the limit of infinite c quark mass this decay of a = l"*", j = 3/2 state should proceed 
via a pure D-wave j^l]. The corresponding decay of its partner, the Dsj(2460)"'", which is 
believed to be a 1"*", j = 1/2 state is energetically forbidden, but if it were allowed it would 
proceed via a pure S wave. Since heavy quark symmetry is not exact, the two 1^ states 
can mix with each other. In particular, an S wave component can appear in the decay 
Dsi{2b3Q)^ — > D*K. Moreover, even if the mixing is small, the S wave component can give 
a sizeable contribution to the width because the D wave contribution is strongly suppressed 
by the small energy release in the Dsi(2536)^ D*K decay. The angular decomposition in 
S and D waves for the analogous decays of the 1"^, j = 3/2 Di{2A20)^'^ mesons to D*^7i^ , 
D*^7i~^ was performed more than 10 years ago by CLEO [5J, but currently no results on the 
D,i(2536)+ exist. 

This study is based on a data sample of 253 fb~^ collected at the T(4S') resonance and 
28 fb"^ at an energy 60 MeV below the resonance with the Belle detector at the KEKB 
asymmetric-energy e~^e~ (3.5 on 8 GeV) collider Q. The Belle detector is a large-solid- 
angle magnetic spectrometer that consists of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a 
barrel- like arrangement of time-of-flight scintillation counters (TOF), and an electromagnetic 
calorimeter comprised of CsI(Tl) crystals (ECL) located inside a super-conducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux-return located outside of the coil is 
instrumented to detect mesons and to identify muons (KLM). The detector is described 
in detail elsewhere 0. Two inner detector configurations were used. A 2.0 cm beampipe 
and a 3-layer silicon vertex detector was used for the first sample of 155 fb~^, while a 1.5 cm 
beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were used to record 
the remaining 126 fb~^ [s^. 

In Monte Carlo, Z^si(2536)^ from e'^e~ annihilation, particle decays and the detailed 
detector response is simulated using the PYTHIA, EvtGen and GEANT packages respectively. 
The and decay modes used in reconstruction are generated with their resonant 
substructures taken from the PDG [iq| but neglecting any interference effects. The width of 
the Dsi{2536)~^ resonance in the simulation is set to zero. Only the D wave matrix element is 
used for the Dsi(2536)^ — > D*~^Kg decay. As is shown below no clear resonant substructure 
is visible in the decay Dsi(2536)^ D^n^K^. Therefore, it is simulated as a three-body 
phase space decay. 
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D,i(2536)+ ^ D+TX-K+ DECAY AND CALCULATION OF A+""^nr+ ^o!^ • 

vr^ and K"^ candidates are required to originate from the vicinity of the event depen- 
dent interaction point. To identify kaons, the dE/dx, time of flight and Cerenkov hght 
yield information for each track are combined to form kaon Ck and pion C.,^ likelihoods 
and the requirement LkI{C.k + l^i^ > 0.1 is imposed. Kg candidates are reconstructed 
via the tt'^tt~ decay channel. and mesons are reconstructed using K~it~^, K^tt+tt", 
K~n'^7i^7i^ and Kg7i~^, K~'K^Tt^ decay modes, respectively. All combinations with masses 
within ±20 MeV/c^ of the nominal D mass are selected and then a mass and vertex con- 
strained fit is applied. 

Candidate D*^^^ are reconstructed using the D^n^ mode. The slow 7r+ momentum 
resolution suffers from multiple scattering. It is improved by a track refit procedure in 
which the vr"*" origin point is constrained by the intersection of the D momentum and the 
known region of e~^e~ interactions. The same procedure is applied to the slow vr" and 
from the Dsi(2536)+ — D^tt~K'^ decay. The M£)0^+ mass is required to be within 
±1.5 MeV/c^ around the D*~^ nominal value. The D*~^ mass constraint fit is not applied. 
Instead, the mass difference M^jo^+^o — Mj;)0^+ is used for Dsi(2536)+ where the errors in 
D*~^ momentum essentially cancel out. 

It is known that the momentum spectrum of the excited charm resonances from contin- 
uum e+e~ annihilation is hard. In addition, due to the strong magnetic field in BELLE 
the reconstruction efficiency for slow n"^ and mesons is larger for higher Dsi{2536)~^ 
momenta. Therefore, to reduce the combinatorial background, it is required that xp > 0.8, 
where the scaled momentum xp is defined as xp = p*/Pmax- Here p* is the momentum of 
the Dsi(2536)"*' candidate in the e~^e~ center of mass frame, p^ax is the momentum which 
the candidate would have if it carried all of the beam energy -Ebeam same frame: 



Pmax= V^Seam-M2. 

The mass M£)+^-k+ and the mass difference {Mj^o^+^o — Mp,o^+) + M]^^^ for all accepted 
combinations are plotted in Fig. [2 Here and in the following the PDG superscript denotes 
the nominal mass value from [l^. A clear peak for the new decay channel Dsi{2536)~^ 
D^Tx~K^ is visible in Fig. H The mass spectrum of the wrong sign combinations D~^7[^K~ 
shown by the hatched histogram has no enhancement in the Dsi(2536)"'" region. 

To calculate the number of Dsi(2536)^ decays, the distributions in Fig. Q are fit to 
the sum of two Gaussians. Their central values are required to be equal. To ensure that 
the second Gaussian is always wider than the first one, its width is chosen to be of the 



form (72 = ycrf + Aa'^. The position of the peak, cti. Act, the fraction of events in the 
first Gaussian and the total number of events in two Gaussians are allowed to vary in the 
fit. The background for the three-body D~^n~K~^ mode is parameterized by the second 
order polynomial muhiplied by the function {M - M%ll'l'°j^if , where M%ll'!^°jli = Mj^'^^ + 
M^P^+M^+^ . For the two-body D*^Kg mode the background parametrization is chosen to 
be of the form — M^J+^o times a first order polynomial, where M^J"+^o°''^ = Mj^B^ + 
M^B'~^ . Table U contains the fit results together with the parameters of the Gaussians 
obtained from Monte Carlo. There is a small fraction of events that contribute two entries 
to the Dsi(2536)"'" signal region in the mass plot. The last row in the Table El shows the 
excess of such events in comparison with the same number averaged over the left and the 
right sideband. The signal and the sidebands are defined as |AMr)+ | < 5 MeV/c^, 10 

si 
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FIG. 1: i:>si(2536)+ mass spectra for D+tt~K+ (top) and D*+K'^ (bottom) decay modes. The 
latter is calculated using the mass difference {M^Ot^+xo — Mjjo^+) + M^^^ . The hatched histogram 
in the top plot shows the corresponding spectrum of wrong sign D~^Tr~^K~ combinations. The fit 
is described in the text. The fit results are listed in Tabled 

MeV/c^ < |AMr,+ I < 20 MeV/c^, respectively, where AMr-)+ is measured relative to the 

si si 

peak position obtained from the fit to Fig. [TJ 



TABLE I: Fit results of the spectra shown in Fig. ^ number of events in two Gaussians, fraction 
of events in the narrower 1st Gaussian, ai, Aa and Mj)^^ — M^^'^. The last three values are 
given in MeV/c^. Note that the value of MJ^^^'^ is known with errors of ±0.34 ± 0.5 MeV/c^. The 

sigma of the 2nd Gaussian is cj2 = \J a\ + Acr^. The 3rd and 5th columns contain the corresponding 
Gaussian parameters obtained in Monte Carlo. The number of double counted events after sideband 
subtraction is given in the last row. 





{D+^-K+)Data 


{D+7i-K+)mc 


{D*+Kl)Data 


{D*+KI)mc 


N events 


802 ± 56 




3474 ± 64 




1st G. fraction 


0.58 ±0.07 


0.769 ±0.014 


0.60 ±0.04 


0.862 ±0.013 


0"! 


0.79 ±0.07 


0.478 ±0.011 


0.98 ±0.04 


0.746 ± 0.014 


Act 


3.1 ±0.6 


1.85 ±0.11 


2.46 ±0.15 


2.49 ±0.18 




-0.51 ±0.06 


0.014 ± 0.008 


-0.42 ±0.03 


0.023 ±0.012 


double counting 


28 - f = 18 




105 - ^ = 96 





To cross-check the results, the mass spectrum is plotted in Fig. |21for the i5si(2536)"'" 
signal and sidebands. The latter is normalized to the area of the signal interval. The 
sideband subtracted plot shown in the bottom of Fig. |21is fit to a double Gaussian as above 
and a constant background. The resulting yield 739 ±51 is consistent with the yield 679 ±48 
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obtained from the fit of the Dsi{2536)~^ mass spectrum. The constant background level is 
found to be —0.4 ± 0.6, which is consistent with zero. The enhancement in the mass 
region observed in the Z)si(2536)+ sidebands (top plot of Fig. |21) is due to combinations of 
a real with a random ti" pair in the event. 



rn Dsl signal 
side bands 



1860 1880 
M(D"'7r"K"') 




1920 



1920 
MeV/c^ 



|AM£,+^-^+| < 5 MeV/c^, open 
normalized to 



FIG. 2: mass spectrum for the Dsi(2536)^ signal band 
histogram in the top plot) and the sidebands (10 < |AM£)+^-/^+| < 20 MeV/c 
the signal interval, hatched histogram). Here AM£)+^-x+ = M£,+^-x+ 
Dsl (2536)"'' peak position in the top plot of Fig. ^ The bottom plot shows the difference. The 
solid curve shows the results of the fit described in the text. 



M9, is the 



Due to the low momenta of the final state particles, the Dsi(2536)"'" reconstruction ef- 
ficiency strongly depends on its momentum, which is found to be harder in data than in 
Monte Carlo. Therefore the ratio of Z)si(2536)"'" branching fractions is calculated using the 
following formula: 

ig(£',i(2536)+ ^ D+TI-K+) _ 
i3(D,i(2536)+ ^ D*+K^) ~ 
-1 

B{D*+^ L)07r+)g(D°)g(ir°-> ^ tt+tt-) 
B{D+) ■ 

Here, e5+^-j^+ and e^.+^o are the reconstruction efficiencies in individual momentum bins 

of Dsi(2536)"'", N^,+j^o are the number of Dsi{2536)~^ D*^Kg decays observed in a given 

momentum bin, and the sum runs over the momentum bins with xp > 0.8. Np)^^^ d*+kI 

and ^£)^''x"£).+;^o are the total number of decays obtained by performing the fit to the 

(2536)"*" mass spectra and the number of double counted events, respectively (see Ta- 
ble H}. B{D^, D^) is the sum of branching fractions of and modes used in the 



„dbl cnf] 

DnK - rij^^j^ J 
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reconstruction For modes with K^^s in the final state the branching fractions are in 
addition multiplied by B{K° K'^ vr+vr ) = 2 ■ B{K^ t^^t^ )• The dependence of the 
efficiency on the initial polarization of the Dsi(2536)^ is checked and found to be negligible. 
The efficiency does not depend on the D*^ helicity angle in the D*^K^ decay and thus on 
the proportions of S and D waves. For the D^7i^K~^ mode the efficiency is independent of 
the D~^TT' , K^TT^ masses and the D^n^ helicity angle. The ratio of branching fractions is 
found to be 

The systematic error receives contribution from different sources listed in Table ITTl A pos- 
sible difference between the data and Mone Carlo in evaluation of the tracking efficiency 
was estimated using partially reconstructed D*~^ decays. The tracking errors due to slow 
K~^, TT^ and pions from Kg are added linearly. The uncertainty in the kaon particle identi- 
fication is estimated using D*~^ decays. Uncertainty in the ratio of and efficiencies 
was determined by a comparison of different decay modes used in the reconstruction. The 
largest contribution to the systematic uncertainty arises due to the fitting model and was 
evaluated by comparing the fit results using different binnings of the D^tt~K~^ and D*~^Kg 
mass spectra. The largest discrepancy (in the D'^n'K'^ decay mode) is taken as the sys- 
tematic error. Finally, the contribution due to the assumption that the efficiency does not 
depend on the decay angles of the Dsi{2536)'^ and its initial polarization, Dtt helicity angles 
and M{D'^n~) is evaluated by comparing the yields of events using either an average or 
differential efficiency in the specified variables. The total systematic error is found to be 
12% (Table HH). 

TABLE II: Systematic uncertainties for ^j^£)^^^2536)+-^D*+K^^) ^ • 

Source Uncertainty, % 

Slow vr^ tracking efficiency 1.5 
Slow tracking efficiency 1 
Slow Kg tracking efficiency 5 

Slow K^ particle identification 1.2 

Ratio of and efficiencies 3.5 

Fit to Af(L>+7r"K+) distribution 8.5 
Efficiency independence on Dsi{2536)^ polarization, 

Dtt helicity angle and M{D~^7r~) 1.1 

Total 12.0 



The D^n^ and K^n^ mass distributions for the Z^si(2536)^ D^ti^K^ decay are 
shown in Fig. |21 The Dsi{2536)~^ signal yield is obtained from fits to the D^tt^K^ mass 
distribution in bins of D^n^ and K^tt~ mass. All Gaussian parameters except the total 
number of events are fixed in the fit to the values obtained from Fig. ^ The position of 
the threshold used for the background description depends on the chosen bin. The areas 
of the plots have been normalized to unity. The plots are not efficiency corrected since the 
differential efficiency does not depend on D^n^ or K^7r~ masses to within the errors of 
Monte Carlo statistics. No dominant resonant substructure is visible in Fig. El 
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FIG. 3: Normalized mass spectra of D~^it~ (top) and K'^tt~ (bottom) pairs from Dsi{2536)~^ 
D~^Tr~K~^ decay obtained from fits to the D~^ir~K^ mass distributions in different D~^Tr~ or K~^Tr~ 
mass bins. The dashed histograms show the corresponding distributions for the Dsi(2536)+ phase 
space decay in Monte Carlo. 

ANGULAR ANALYSIS OF Dsi(2536)+ ^ D*+K'^ DECAY. 

The Dsi{2536)~^ D*~^Kg decay can be described by three angles a, (3 and 7 defined 
as shown in Fig. ^ The angles a and (3 are measured in the Dfi rest frame with respect to 
the direction of the boost needed to go from the e'^e' center of mass frame to the Dfi rest 
frame, a is the angle between the boost direction and the Kg momentum. (3 is the angle 
between the decay plane and the plane formed by the Kg and the boost direction. The third 
angle 7 is defined in the D*~^ rest frame between vr"*" and Kg. 




K 



FIG. 4: Definitions of the angles a, (3 and 7. The first two are measured in the Dsi(2536)^ rest 
frame, the last one - in the D*^ frame. "Boost" changes the e^e~ center of mass system to the 
L>si(2536)+ frame. 

The measured cos a, (3 and cos 7 angular distributions are shown in Fig. El They represent 
the signal yield obtained from fits to the D*^K^ mass spectra in bins of individual angular 
variable. All Gaussian parameters except the normalization are fixed in the fit to the values 
obtained from the overall spectrum in Fig. ^ The first distribution is efficiency corrected 
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since there is a slight hnear efficiency dependence on cos a: the relative difference between the 
values at —1 and at +1 is about 12%. Since the differential distribution in cos a includes only 
even powers of cos a, such a linear dependence does not induce any biases when integrating 
over the whole allowed interval [—1, +1] and thus no dependence of the integrated efficiency 
on the initial polarization of Z)si(2536)"'". No significant efficiency dependence on (3 and cos 7 
is observed and therefore the corresponding two spectra in Fig. El are not efficiency corrected. 
One can see that the first two distributions are not fiat. This means that Dsi(2536)"'" is 
produced polarized and that it does not decay in a pure S wave. The last distribution 
is more important. For a pure S or D wave decay it should either be fiat or have the 
form (1 + 3cos^7), respectively. In the general case of interference between S and D wave 
amplitudes it becomes a linear combination of sin^ 7 and cos^ 7: 

l^T^ = ^{^+ (1 - R){^-^^^) + \/2i?(l-^)cos0(l - 3cos^)}, 

where R = Ts/ (Ts + Td), Ts,d are the S, D wave partial widths respectively, is the relative 
phase between the two amplitudes. This is similar to the case of Djj°(2420)*^ —>■ D*^7i~ and 
D*(2420)+ D*°7i+ decays studied by CLEO 5]. 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 
cos Y 



FIG. 5: Angular distributions for the -Dsi(2536)^ D*~^Kg channel. The definitions of angles a, 
f3 and 7 are given in the text and in Fig. |1J Dsi{2536)^ reconstruction efficiency depends slightly 
on cos a, therefore the upper left plot is efficiency corrected. A fit to the cos 7 distribution is 
described in the text. 



Fitting the cos 7 distribution in Fig. El to the form 1 + 74cos^7 yields A = —0.70 ± 0.03 
and a P^^^ degree of freedom of 1.39. It was checked that the A parameter obtained on 
the subsamples defined by the cuts cos a < 0, cos a > 0, P < 71 /2, [3 > 7r/2, xp < 0.93 or 
xp > 0.93 agree with each other within statistical accuracy. Knowledge of A constrains the 
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contribution of S wave to the width and the relative phase 0: 

3-A _ p 

2^2R{1 - R) 

The allowed range |cos</)| < 1 is shown in Fig. El Two lines in this plot bound the re- 
gion which corresponds to ±1ct deviation in A. Regardless of the value of 0, the S wave 
contribution is limited from below and from above by the values corresponding to = 0. 
Conservatively taking the value of A to be 2(T below the central value, the limits obtained 
are: 

0.277 <R< 0.955. 
The corresponding limit for is: |0| < 42°. 




'o 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
R 



FIG. 6: Plot of cosine of the relative phase of S and D wave amplitudes in the -0^1(2536)''' D*~^K^ 
decay versus R = r5/(r5' + To)- The two curves bound the region that corresponds to a ztlo" 
deviation in the measured parameter A. 



In conclusion, a new decay channel Dsi(2536)^ D^tt^K^ is observed. The D+7r~ 
pair is the only Dn combination that cannot come from a real D* resonance. It can be 
produced in Dsi(2536)"'" two-body decays only through the virtual resonances D*^, broad 
Dq^ or D2°(2460). In addition, the D~^7r~K~^ final state can be formed by two-body decays 
to a and a virtual K*^ or higher K* resonance. No clear resonant substructure is found in 
the D~^7c~K~^ system. The ratio of branching fractions ^^[^ f(2536|+-fz? *+ k'^) ^ measured to 
be (2.8±0.2±0.4)%. An angular analysis of the normalization decay Dsi(2536)+ D*~^Kg 
is also performed. The Dsi{2536)~^ may mix with another = l"*", j = 1/2 state, which is 
presumably the recently discovered Dsj{2460)~^ meson and can decay in an S wave. Since 
the energy release in this reaction is small, the D wave is suppressed and the S wave can 
give a sizeable contribution to the total width even if the mixing is small. The measured 
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1 — (0.70 ± 0.03) cos^7 D*^ helicity angular distribution constrains the relative fraction of 
the S wave component to the range 0.277 < R < 0.955, independent of the phase 0. 
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